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Abstract 

We calculate the linear and circular polarization of gluons produced in con- 
junction with massive quarks in the annihilation process e^e~ — > qqG. The 
linear polarization is calculated in the hadron event plane as well as in the 
gluon-beam plane. Beam polarization and polar orientation effects are in- 
cluded in our discussion. For typical top pair production energies at the 
Next-Linear-Collider (NLC) the degree of linear polarization in the hadron 
event plane remains close to its soft gluon value of 100% over most of the en- 
ergy spectrum of the gluon. The linear polarization in the gluon-beam plane 
is generally smaller but peaks toward the hard end of the gluon spectrum. 
The dependence of the linear polarization on beam polarization and on the 
polar orientation of the gluon is small. The circular polarization is largest for 
maximal gluon energies and shows a strong dependence on the longitudinal 
beam polarization. The longitudinal polarization of the beam may therefore 
be used to tune the circular polarization of the gluon. The massive quark 
results are compared with the corresponding results for the massless quark 
case. 



^on leave absence from GIF, Colombia 



1 Introduction 



The polarization of gluons in e^e~ annihilation |jT|, |], in deep inelastic scattering ||^ 
and in quarkonium decays Q has been studied in a series of papers dating back 
to the early 1980's. Several proposals have been put forward to measure the po- 
larization of the gluon among which is the proposal to measure azimuthal angular 
correlation effects in the splitting process of a polarized gluon into a pair of gluons 
or quarks f^. Latter proposal has led to a beautiful confirmation of the presence of 
the three-gluon vertex using e^e~ data (see also 

The earlier calculation of the gluon's polarization in e^e~ annihilations had been 
done for massless fermions which was quite sufficient for the purposes of that pe- 
riod. In the meantime the situation has changed in as much as the heavy top quark 
has been discovered at Fermilab in 1995 whose production properties in e^e~ an- 
nihilations will be studied in the proposed Next-Linear-Collider (NLC). In its first 
stage typical running energies of the NLC would extend from tt-threshold at about 
350 GeV to maximal energies of about 550 GeV. It is quite clear that top mass 
effects cannot be neglected in this energy range even at the highest cm. energies. 
It is therefore timely to redo the calculations of []1|, for heavy quarks and to inves- 
tigate the influence of heavy quark mass effects on the polarization observables of 
the gluon. A first step in this direction was taken by us in [^] where we determined 
the linear polarization of the gluon in the process e^e~ ttG taking the hadron 
event plane (for short: event plane) as a reference plane. In |p] we did not take 
into account azimuthal and polar orientation effects of the event plane relative to 
the beam. This is done by an appropriate polar and azimuthal averaging process. 
Beam polarization effects were only commented on in passing in p. 

In this paper we extend the analysis of [§| in several directions. We include beam- 
event orientation and beam polarization effects in our discussion. We also compute 
the circular polarization of the gluon induced by the parity-odd component of the 
hadron tensor and/or by longitudinal beam polarization effects. Finally we compute 
the linear polarization of the gluon in the gluon-beam plane which is obtained after 
an appropriate azimuthal averaging process. 

Our paper is structured as follows. Sec. 2 describes the general formalism of 
massive quark plus gluon production in e^e~ annihilations including beam and 
gluon polarization effects. Our production cross section is written in terms of three 
modular building blocks. The first building block defines the orientation of the 
lepton beam relative to the hadron plane and the dependence on the polarization 
parameters of the beam. The second building block specifies the electro-weak model 
dependence and the third building block specifies the QCD dynamics in terms of a 
set of polarized and unpolarized structure functions. In Sec. 3 we list our results for 
the twice-differential 0{as) polarized and unpolarized structure functions as well as 
closed form expressions for their once- and twice-integrated forms. Sec. 4 contains 
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our numerical results. We provide plots of the energy dependence of the linear 
and circular polarization of the gluon, their polar angle dependence and their beam 
polarization dependence. In Sec. 5 we discuss the linear polarization of the gluon 
in the gluon-beam plane which can be obtained via an azimuthal rotation from the 
event plane. Sec. 6 contains our summary and our conclusions. 

Before we get to the main topic of this paper we want to briefly present some 
numerical evidence for what is referred to as the "dead-cone" effect which occurs 
when soft gluons are radiated off heavy quarks. In the heavy quark case there is 
a depletion of gluons close to the direction of the heavy quark and the antiquark. 
This is quite different from the case of gluons radiated off light quarks where the 
production peaks towards the collinear limit. Large angle emission is a welcome 
effect in e~^e~ annihilation since gluons emitted at large angles from the heavy 
quark and antiquark directions are easier to reconstruct. If one is far enough away 
from the production threshold with sufficient velocities of the heavy quark and 
antiquark the decay products of the heavy quark and antiquark would tend to go 
along the original production direction and would thus stay away from the large- 
angle soft gluons radiated from the heavy quark and antiquark. For example, at 
= 500 GeV and ^/q^ = 1000 GeV one has a quark (or antiquark) velocity of 
V = 0.71 and v = 0.94, respectively. 

In order to display the "dead-cone" effect we have plotted the quark/gluon open- 
ing angle distribution at a cm. energy of = 500 GeV for four different soft gluon 
energies in Fig. 1(a). The quark gluon opening angle ^13 is given by 

+ {2-x)w 

cos alls = — , (1) 

At these low energies there is in fact a depletion of gluons radiated at small opening 
angles. However, the distributions rise uniformly to their maximum values at 90° 
and do not show a peak at around 6*13 = 2m/y/q^ = 40.10° as predicted in 



In Fig. 1(b) we show the opening angle distribution for gluons at a fixed energy of 
E = 6 GeV for three different values of the cm. energy. Again there is a depletion 
of gluons radiated at small opening angles for all three cm. energies. The \/q^ = 
1000 GeV opening angle distribution does show a peak at around di^ = 40° which 
is somewhat displaced from the peak position 6'i3 = 20.05° as calculated from the 
above relation. The two distributions at lower cm. energies do not show any peak 
structure but rise uniformly to their maximum values at 90°. 

We mention that some indirect evidence for the "dead-cone" effect was presented 
in |Tl| where we calculated the mean transverse momentum of the top with regard 
to the antitop direction in e'^e~ annihilations. The mean transverse momentum of 
the top was found to be only slightly below the mean gluon energy in the whole 
range from tt-threshold to 1000 GeV. Since the transverse momentum of the top 
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has to be balanced by the transverse momentum of the gluon the near equahty of 
the two means imphes large mean opening angles for the gluon. 

2 General formalism 

As is usual we shall represent the two-by-two differential density matrix dux = dcxx^x'^ 
of the gluon with gluon helicities Xq = ±1 in terms of its components along the unit 
matrix and the three Pauli matrices. Accordingly one has 

dsT = ^{dal + da^'sT^ + da^sTy + rfa^dr^), (2) 

where da is the unpolarized differential rate and da = {da^, da^, da^) are the three 
components of the (unnormalized) differential Stokes vector. In this paper the 
{x, y, 2;)-components of the Stokes vector will mostly be referred to the event plane 
such that the gluon points in the ^r-direction and the (ggG)-plane defines the (x, z)- 
plane (for short: event plane). An exception is Sec. 5 where the linear polarization 
of the gluon is calculated in the gluon-beam plane. 

Specifying to e~^e~ — > q{pi)q{p2)G{p3) we perform an azimuthal averaging over 
the relative beam-event orientation. After azimuthal averaging the ^/-component of 
the Stokes vector da^ drops out |l|, One retains only the x- and z-components 
of the Stokes vector which are referred to as the gluon's linear polarization in the 
event plane and the circular polarization of the gluon, respectively. The differential 
unpolarized and polarized rates, differential with regard to the polar beam-event 
orientation and the two energy- type variables x = 2p3 ■ g/g^ and w = 2(pi — ^2) ■ g/g^ 
(with q = Pi + P2 + Pd) are then given by 

d cos 9 dx dw 



da^ 
d cos 9 dx dw 

The notation da^^^ stands for either da or da^, and the same for rfcr^^^^ (the indices 
i = 1,2,4 and a = U,L,F are explained later on). The notation closely follows 
the one used in 0. Noteworthy is the absence of a longitudinal contribution to the 

^The linear gluon-beam plane polarization observable of the gluon discussed in Sec. 5 involves 
a different azimuthal averaging process. For this application one needs to retain one particular 
y-component of the Stokes vector in the event plane. 
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^ [1 + cos' 9) i gj""^ 



l{x) 



dx dw 



+ 912 



da 



2{x) 



dx dw 



(3) 



4 x^^^dxdw 



rfafM 3 , daf^ 
dx dw / 4 dx dw 



|(1 + cos^»)9,,i^ + 5 „os» + . 

8 dx dw 4 \ dx dw dx dw I 



(4) 



circular polarization with an angular sin^ ^-dependence. This is a tree- level effect 
due to the CP-evenness of the Standard Model interactions. 

We have written the electro-weak cross section in modular form in terms of 
three building blocks. The first building block determines the angular beam-event 
dependence. For the case at hand one remains with a polar angle dependence after 
azimuthal integration, where 9 is the polar angle between the gluon and the elec- 
tron beam. The second building block specifies the clcctro-weak model dependence 
through the parameters gij = 1, . . . , 4). They are given by 

911 = Q) -2Q}VeVfRexz + + a^^iv) + a))\xz?, 

912 = Q) -2QfVeVfRexz + {vl + al){v) - a))\xz\^, 

9u = 2QfVeafRexz - ^{vl + al)vfaf\xz\^, (5) 

941 = ^Qja^VfRexz -2vea^{v) + a))\xz?, 

942 = 2Qfa,A}fRcxz -2veae{v) - a))\xz\^, 

944 = -2QfaeafRexz + 4:VeaeVfaf\xzf, 

where, in the Standard Model, Xz(g^) = 9Mlq^/{q^ + iMzTz), with Mz and 

Tz the mass and width of the Z° and g = GpiSv^Tra)-^ ^ 4.49-10-^ GeVl Qf are 
the charges of the final state quarks; fg and Cg, f/ and af are the electro- weak vector 
and axial vector coupling constants. For example, in the Weinberg-Salam model, 
one has Vg — —1 + 4sin^^i^, — —1 for leptons, Vf — 1 — Isin^^^^, af — 1 for 
up-type quarks [Qf — |), and Vf = —1 + |sin^^iy, a/ = —1 for down-type quarks 
[Qj = -|). In this paper we use Standard Model couplings with sin^ 6y/ = 0.226. 

The third building block, finally, is given by the hadron dynamics, i.e. by the 
current-induced production of a heavy quark pair with subsequent gluon emission. 
The QCD dynamics is encoded in the hadronic rate functions da'^^^y dx dw and 
dal^/dxdw with specific components denoted by the indices a and i, and the polar- 
ization indices x and z. The hadronic rate functions can in turn be related to the 
polarized and unpolarized hadronic structure functions H^^^^ and according to 



da'S^^ . dai 



„2 



The different hadronic structure functions are specific components of the hadronic 
three-body tensor Hj^i, induced by the product of vector current (V) and axial vector 
current (A) interactions. In the case of massive quark pair production one has 
four different independent current products which are denoted by the upper index 
i — 1, . . . ,A, where 

2^ 2^ ' 
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We have temporarilly dropped all further indices on the hadron tensor in Eq. (|^). 

We mention that, because of CP- invariance and because we are calculating at 
tree level, the linear combination does not contribute to the cross section and 
the spin observables considered in this paper. The lower index a = U, L, F specifies 
the relevant density matrix element of the gauge boson that determines the polar 
angle cos6'-distribution. The three elements of the density matrix are referred to as 
the "unpolarized transverse" component Hu = H^+ + H = H11 + H22, the "longi- 
tudinal" component = Hqq and the "forward-backward asymmetric" component 

Hp = — H = —i{Hi2 — H21) where we have listed both the spherical and 

Cartesian components of the relevant density matrix elements. The three (?7, L, F) 
components can be obtained from the full hadron tensor H^^, by covariant contrac- 
tion as discussed in 0. The unpolarized hadronic tensor and the polarized 
hadronic tensors H^^ and H^^ are defined in exact analogy to Eq. (^). 

Eqs. and give the differential cross section for unpolarized beams. The 
case of longitudinally polarized beams can easily be included. For the unpolarized 
and linearly polarized rates da]^^^ one has to effect the replacement 

gii ^ (1 - h-h+)gii + {h- - h+)gii {i = 1, 2) 

944 {h- -h+)gi4 + il-h-h+)g44 (8) 

where /i+ and /i^ are the degrees of the longitudinal (or helicity) polarization of the 
positron and electron beam, respectively. For the circularly polarized rates da^^ one 
has the replacements 

5^14 ^ (1 - h~h^)gi4 + {h' - h^)g44 

g4i ^ {h- ~h+)gu + {l-h-h-^)g4^ (^ = 1,2). (9) 

Transverse beam polarization effects can also be easily included in the framework of 
our formalism (see e.g. P) but will not be further discussed in this paper. 

After having described the general formalism we now go on to present explicit 
expressions for the polarized and unpolarized hadronic structure functions for the 
annihilation process e^e" qqG in the following section. 



3 Polarized and unpolarized structure functions 

The various components of the hadronic tensor can be easily calculated from the 
relevant tree level Feynman diagrams and are given by 

2 1 2 

H}j(x,w) = -16(- -2 + x)- -6At+(x,w) + 32(- -2 + x)tUx,w), 

X ^ X X 

Hl{x,w) = -32^U{x,w) + 16e(^ - 2y4x,w), 
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HUx,w) = 32f--l)- + 16^t+(x,w;)-8^f - 2 - x)tUx,w), (10) 

2 c 

H^{x,w) = -16^t+(x, w) + 8n ^ - 2 - xH^ix, w), 

Hi(x,w) = 64t(x,w)-32(--2 + xY-{x,w), 

^x ' 

E]^{x,vS) =- 32(--l)--64t+(x,w) + 64(--l)4(x,w), 

HIj%x,w) = -32^t+(x,w) + 16^(^^-2)4(a;,u7), 

H}''(x,w) = 32(- -l)- + lQ^t+(x,w) -8^(^^ -2)tUx,w), (11) 

i/f (x, w) = -IG^Uix, w) + 8e(^ - ^yHx, w), 

Hp(x,w) = 64t_(x,M;) -64f- -l)ti(x,?x;), 

^x ' 

Hl;{x,w) = 64xt„(x,w) -32(2-x)ti(x,w), 
Hl^{x,w) = 0, 

H},^{x,w) = -16{2-x)--64:Xt+{x,w) + 32{2-x)tUx,w), (12) 

a; 

Hf{x,w) = 
with ^ = 4m^/g^, where we have used the abbreviations 

t±(x, w) := ^ ( ^ ± ^ ) , ti(x, ti;) := ^ ± ^— ) . 

4 Y(x — w)^ (x + wj^y 2 Vx — w x + wj 

(13) 

We mention that the vanishing of the hadron tensor components Hf^^{x,w) and 
Hp^{x,w) is a tree-level effect and is due to the CP-invariance of the underlying 
Standard Model interaction. 

Note that under quark-antiquark exchange (charge conjugation) t± ±t± and 
— i> ±t^. This implies that the charge conjugation odd contributions t_ and ti_ 
vanish if one does not discriminate between quarks and antiquarks. In the following 
we wish to distinguish between the cases where the quark flavours are identified or 
not identified and we will refer to these two cases as the flavour tag and flavour 
no-tag cases, respectively. 

Let us briefly pause to discuss the mass zero limit of the hadronic tensor com- 
ponents. The mass zero limit can easily be taken by setting ^ = in Eqs. ([T0|-p!2D. 
From the terms that remain after taking the ^ — > limit it is only the terms propor- 
tional to t^^{x,w) that are important, because they are mass singular in the ^ — > 
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limit. In fact, when one performs the w-integration including flavour-tagging, the 
mass singluar functions ti_{x,w) integrate to an x-independent logarithmic factor 
and a finite (x-dependent) term, i.e. / tj.{x,w)dw — >• — In^ + c± (see Eq. (p!3|)). 
Keeping only the dominant logarithmic term one has in the ^ — >^ limit 

2 2 
Hhix) -32(- - 2 + x) In^, H^{x) 32(- - 2 + x) In^, 



HlTix) ^ -64(i - l) Ine, Hp{x) ^ 64(^ - l) In^, (14) 
Hl;{x) 32(2 - x) In^, i^)?(x) ^ -32(2 - x) In^. 



The remaining terms in Eqs. (|T0|-|l2|) are subdominant. The expressions in Eqs. (|T^ 
determine the energy dependence of the linear and circular polarization of the gluon 
in the mass zero limit. They will be used later on to compare with the corresponding 
expressions in the massive case. 

Returning to the case of massive quarks, the integration over the phase-space 
parameter w can be easily done. However, as mentioned earlier on, C-odd observ- 
ables average out to zero in this integration if one does not employ flavour-tagging. 
For example, the contributions to the circular polarization vanish when we integrate 
over the whole range of w since the integration variable is antisymmetric in the 
quark and antiquark variables. Nonvanishing values of the circular polarization are 
obtained only if one applies flavour-tagging, i.e. one has to take care to distinguish 
between quark and antiquark energies in the integration. Thus, for a flxed value 
of the gluon's energy, one flrst integrates over positive values of the w-variable and 
then subtracts the integral over negative w- values. Or, equivalently, one just takes 
twice the value of the circular polarization calculated for positive w. 

Let us list the few basic w-integrals that are needed in the calculation. The 
integration has to be done between the symmetric phase space boundaries — w+ and 
w+ where 

w+{x) = xJ^—^-^. (15) 



One has 



t+{x,w)dw = ( l)— w+(x). 



ui+(a;) ^X ' X 



t^{x,w)dw = In =: t+ x), 16 

■w+{x) - X - ^yl - X - t, J 

r+w+ix) 1 — ^ 

2 / t-{x, w)dw = 1, 

JO X 

2/ tt{x,w)dw = In p— ^ =: ti(x). 
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where we have employed flavour-tagging for the last two C-odd integrals as discussed 
before. Using these basic integrals we obtain the once-integrated structure functions 



if^(x), H^^{x) and H^^{x). One has 





= -32(-- 

^x 




^x) 


= -32e(-- 

^x 




^x) 




Hh 




^x 


Hp 


' x^ 


= 64(1 -« 

^ X 




» 


/I 

= -128(-- 

^x 


HlTi 


[x) 


= -32e(-- 

^x 


Hpi 


[x) 


= 16(4 + f) 


Hpi 


[x) 


^x 


Hp{ 


[x) 


= 64(1 

^ X 


Hpi 


[x) 


= 64(1 -^ 


Hpi 


[x) 


= 0, 


Hp{ 


[x) 


= -32(4 - , 


Hp{ 


» 


= 0. 



1 2 

I + x) -w+(x) + 32(- - 2 + x)tUx), 

l)^u;+(x) + 16e(^-2)4(a;), 
I - l)^t^+(a:) - 8e(^ - 2 - x)4(a:), (17) 

l) ^ti;+(x) + 8^ (^-^ - 2 - x)4(a;), 
- l) -32(^ -2 + x)ti(x), 



X ^x 



X ^ X 



-2]t\ix) 



lA-' iXj 

-l)^«;+(x) + 8e(^-2)t^4x), 
-l)-64(i-l)ti(x), 



x), (18) 



1 

X 



(19) 



Our notation is such that H{x) = J H{x,w)dw. The dominant leading logarithmic 
In^ terms in the zero mass limit derive from the logarithmic functions t'^(x) and 
t'_(x) and have been listed before in Eq. (|I^). 

The last step is the integration over the second phase-space parameter x. It is 
clear that we have to introduce a gluon energy cutoff at the soft end of the gluon 
spectrum in order to keep the rate finite. Denoting the cutoff energy by Ec = A-^/g^ 
the integration extends from x = 2X = 2Ec/ to s = 1 — We obtain 

H^ = -128e?(-l) + 128^?(0)-32e?(l) + 64e?,(-l)-64^?^(0) + 32e?,(l), 



9 



= -32e6?(-l) + 32^6^(0) + 16^(2- 06^^-1) -32ee?^(0), 

Hi = 16(4 + oe^(-i)- 16(4 + 06^(0)- 8^(6 -oe^£(-i) + i6ee^£(o) + 8^6;,(i), 

Hi = -lQ^g{-l) + lQ^g{0) + 8a2-OGe{-l)-lQ^gi{0)-8^ge{l), (20) 

H^p = -8(13 - 12^ - - 8Li2(l-0)-481ne- 64(1 -e + lnO In ( — 



V 



H]f = -1286;(-1) + 1286^(0) + 646;,(-l)-646;,(0), 

H^^ = -32e6?(-l) + 32^6^(0) + 16^(2 -O6^K-l)-32ee?K0), 

= 16(4 + 0^(-l)- 16(4 + 06^(0)- 8^(6 -O6?£(-l) + 16e^£(0), (21) 

Hi- = -16e6?(-l) + 16^6^(0) + 8^(2 -O6?£(-l)-16e6^£(0), 

Hp = -64(2(1 -0 -Li2(l -0) -641x1^-64(1 -e + lnOln 



2y 



Hi; = 24(1 -0(3-0 + 48 Ine, 

Hp = 0, (22) 
Hp = -1286?(0) + 966^(1) + 646?,(0)-326?,(l), 
Hp = 
where 

g(m) := I x'^-^w+(x)dx. (23) 

Full analytical results for the integrals have been given before in |^ . Here we employ 
a small A-expansion and retain only the leading In A and the next-to-leading constant 
contributions in the small A-expansion. One has 

.(-1) ^ -ln(i±^)-.ln(^), (24) 

^(0) = (25) 
^(1) = ^(2 + 0^-|(4-01n(^), (26) 
gJm) ■= C ^x"'tUx)dx, (27) 
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= -^ + ^(2-01n([^), (29) 

Giil) = -^(2-0^ + ^(8-8e + 3ain(^). (30) 

We have checked that the above leading and next-to-leading log expansion has suf- 
ficient accuracy for the numerical applications to be discussed in Sec. 4. 

4 Numerical results 

We are now in the position to discuss the various differential distributions of the 
polarization variables with regard to the gluon energy. The QCD dynamics is com- 
pletely specified by the various hadron tensor components assembled in Eqs. (p!0|)- 



(12). In order to exhibit the quark mass and fiavour dependence of the gluon's 
polarization we shall present results for the top, bottom and charm pair production 
cases in some of the figures. We begin our discussion with the linear polarization 
of the gluon in the hadron event plane defined by the plane spanned by the quark, 
antiquark and the gluon. As discussed before the linear polarization of the gluon in 
the event plane is determined by the normalized x-component of the Stokes vector. 
To start with we integrate out the polar angle dependence in the numerator and de- 
nominator of the ratio of polarized and unpolarized cross sections given in Eq. (^. 
One has 

where da^/dx and da/dx are the twice-integrated polarized and unpolarized differen- 
tial cross-sections appearing in Eq. (^. In Fig. 2(a) we show the energy dependence 
of the linear polarization of the gluon produced in association with charm and top 
quark pairs at a cm. energy of \/q^ = 500 GeV. We have chosen to rescale the 
scaled gluon energy x in this plot to its maximal value given by x^ax = 1 — ^■ 
Fig. 2(a) shows that the linear polarization remains close to its maximal value of 
100% at the soft gluon point over a large portion of gluon phase-space. For the top 
quark case the polarization remains above the value in the charm quark case when 
plotted against x/xmaxQ In Fig- 2(b) we show the differential cross section which 
enters in the normalization of the polarization. The differential cross section for glu- 
ons produced in conjunction with top quark pairs can be seen to be approximately 
one order lower than for gluons produced in conjunction with charm quark pairs. 

■^When the linear polarization is compared at equal gluon energies, this relation is reversed. 
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In Fig. 3(a) we display the linear polarization of the gluon for three different 
values of the cm. energy plotted against the scaled energy variable x/xmax- As the 
lowest cm. energy we choose \f(^ = 370 GeV which is far enough above tt-threshold 
for a perturbative treatment to be valid. The linear polarization remains quite close 
to its soft-gluon value of 100% over most of the available phase-space due to the 
fact that there is not much energy available for the gluon so close to threshold. At 
the hard end of the spectrum the linear polarization has to go to zero for the simple 
reason that one can no longer define a hadronic plane in this coUinear configuration. 
As the cm. energy is increased, the linear polarization decreases for a given fixed 
fractional energy value x/xmax- For the largest cm. energy value i/g^ = 1000 GeV 
in Fig. 3(a) the linear polarization is quite close to the zero mass limiting value 
also shown Fig. 3(a). As discussed before, the zero mass result is determined by 
the collinear configurations of the gluon, i.e. by the In^ contributions in Eq. (|14D- 
Using the results of Eq. (0), one obtains^. 



Note that the zero mass linear polarization no longer depends on the electro-weak 
model parameters, on the flavour of the produced quark or on the polarization of 
the beam. The electro-weak model dependence, the flavour and beam polarization 
dependence is the same in the numerator and in the denominator and drops out 
when taking the ratio. It can be checked that the electro-weak model, the flavour 
and beam polarization dependence of the linear polarization is quite weak also in 
the massive case for the cm. energies considered here. By comparing Fig. 2(a) with 
Fig. 3(a) one sees that the linear polarization in the charm quark case is already 
quite close to the zero mass case indicating that the approach to the asymptotic 
value is quite fast. This also explains the closeness of the ^/q^ = 1000 GeV curve to 
the zero mass curve. 

In Fig. 3(b) we show the differential cross section for e"^e^ — tiG for the same 
three cm. energies as in Fig. 3(a). The differential cross section decreases as the cm. 
energy is increased. As expected, the decrease is governed by the usual g^-factor in 
the differential cross section. 

In Fig. 4 we show the circular polarization of the gluon for the same three cm. 
energies as in Fig. 3 including again the zero mass case for up-type quarks. It is 
clear that we are employing flavour tagging since the circular polarization would 
be zero otherwise. The circular polarization is positive in all four cases and peaks 
toward the hard end of the gluon spectrum^. The circular polarization is zero at the 
soft end of the spectrum which is simple to understand since the linear polarization 

^The limiting value of the polarization agrees with the corresponding result Eqs.(4)-(9) in [Q 
(second reference) in the limit xo,S 

'^Note that the positive value holds for Eq > Eq and becomes negative for Eq < Eq. 




(32) 
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already saturates the total polarization bound of 100% at this point. The circular 
polarization remains quite small over most of the range of gluon energies except 
for the zero mass case. The limiting zero mass behaviour is again governed by the 
coUinear configuration determined by the In^ contributions in Eq. ([l^)- One obtains 



P'(x) = - — ^ ^— 33) 

Note that the ^/q^ = 1000 GeV curve in Fig. 4 is still a large way away from the 
zero mass limit indicating that the approach to asymptotia is much slower for the 
circular polarization than what was observed for the linear polarization. 

The x-dependent factor in the mass zero formula Eq. ([33| ) for the circular po- 
larization starts from zero at the soft-gluon point and reaches unity at the hard 
end of the spectrum. Apart from the x-dependent factor the size of the circular 
polarization is governed by the ratio of electroweak factors —Qu/gn- Above the 
top quark threshold the ratio —gn/gu is only mildly g^-dependent but is strongly 
flavour dependent. For example, for = 370, 500 and 1000 GeV one finds 0.273, 
0.265 and 0.257 for the up-type quarks, and 0.678, 0.670 and 0.657 for down-type 
quarks, respectively. From the numerical values of the ratio —gii/gu it is clear that 
one generally obtains much larger values for the circular polarization of down-type 
quarks than for up-type quarks. In Fig. 4(b) we provide a plot of the x-dependence of 
the circular polarization for -y/g^ = 1000 GeV in the zero mass case for up-type and 
down-type quarks. The above numerical results of the ratio gu/gu can be clearly 
identified at the hard end of the spectrum. For the sake of comparison we have 
also plotted the corresponding curve for bottom quarks with = 4.1 GeV. Mass 
effects tend to reduce the circular polarization showing again that the approach to 
the asymptotic formula Eq. ([33| ) is quite slow. 

In Fig. 5(a) we display the beam polarization dependence of the circular polariza- 
tion assuming that the electrons are polarized (the beam polarization dependence 
of the linear polarization is quite weak). The circular polarization of the gluons 
can be seen to be strongly dependent on beam polarization effects. One observes a 
strong polarization transfer effect from the helicity polarization of the electrons to 
the circular polarization of the gluon which can enhance the circular polarization of 
the gluon up to a factor of ten. In Fig. 5(b) we exhibit the beam polarization depen- 
dence of the differential cross section. The differential cross section is enhanced by 
about 50% for negative electron heli cities and reduced by about the same amount 
for positive electron helicities. 

In Fig. 6(a) and Fig. 6(b) we show the cos ^-dependence of the linear and circular 
polarization of the gluon where we remind the reader that 6 is the polar angle of the 
gluon relative to the electron beam. The dependence is shown for a cm. energy of 
^/q^ = 500 GeV and for three different gluon energies. The cos 6'-dependence of the 
linear polarization is quite weak for all three gluon energies but becomes somewhat 
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larger as one moves away from the soft gluon point where the cos ^-distribution 
is flat. In contrast to this, the circular polarization shows a strong asymmetry 
effect and turns from positive to negative values as cos 6' moves from —1 to +1. 
The asymmetry effect becomes larger as the gluon becomes harder. As Fig. 6(c) 
shows, the pronounced asymmetry of the circular polarization results mainly from 
the numerator of the relevant polarization expression. The differential cross section 
shown in Fig. 6(c) does not possess a very pronounced cos 6'-dependence. 

In Fig. 7(a) we show a plot of the average linear polarization of the gluon as a 
function of the cm. energy for three different cutoff values = Aa/^ = 3, 5 and 
10 GeV. Gluon energies of this magnitude are sufficient to make the corresponding 
gluon jets detectable. The average linear polarization of the gluon rises steeply from 
threshold and quickly attains very high values above 95%. The linear polarization 
becomes larger for smaller values of = Aa/^^ and tends to 100% as the cutoff 
parameter tends to zero. The approach to the asymptotic value = 100% is, 
however, rather slow. We also show the average linear polarization of the gluon 
limit when the quark mass goes to zero for a cutoff energy of 3 GeV. The average 
linear polarization is smaller than in the massive case. In Fig. 7(b) we show the 
cutoff dependence of the total cross section as functions of the cm. energy. The 
cross sections rise steeply from threshold and then turn over to their canonical l/q^- 
behaviour. The cutoff dependence of the cross section is still moderate for the three 
chosen cutoff values. 

5 Linear polarization in the gluon-beam plane 

In Sees. 2, 3 and 4 the linear polarization of the gluon has been specified relative to 
the event plane. A measurement of the polarization would require the reconstruction 
of the event plane which may not always be simple. If is much easier to determine 
the gluon-beam plane since the beam coordinates are known a priori. In order to 
determine the (unnormalized) components of the Stokes vector in the gluon-beam 
plane one needs to rotate the event plane Stokes vector around the gluon momentum 
axis by the azimuthal angle x between the the two planes. One obtain^ 

H"" = cos2xH'' - sm2xHy 

Ry' = sin 2x^'= + cos 2xi^^ (34) 

= 

where [H^' , H^', H^') are the components of the Stokes vector in the gluon-beam 
plane. Since one is rotating around the gluon axis, the circular component is not 

^Note that the Stokes vector is not a true vector and thus its transformation behaviour under 
rotations differs from that of a true three-vector. 
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affected by this rotation. We have dropped all further indices on the hadron tensor 
in Eq. (^). 

We shall not write down the full polar cos 9- and azimuthal x-dependence of the 
polarized cross section in the new system but shall rather completely integrate out 
the polar and azimuthal angle dependence. One then obtains 

Hi' = (35) 

where the y'-component H^' drops out after azimuthal averaging. The hadron tensor 
components "T" and "4" on the right hand side are still specified in the event system 
and are given by (see 0) 

H^T = \{Hl^ + H%) = ^{-H^,, + H^,) 

m = -\{Hl_-Wi^) = -\{Hy, + Hy,). (36) 

Again these components can be calculated from the relevant tree-level Feynman 
diagrams. One finds 

X 

H^'''{x,w) = 0. (37) 

Since we are interested in the energy dependence of the linear polarization, we 
also give their once-integrated forms. One has 

H'/ix) = 8(^-±^-A^-A^^tiix) 

\ X J X X 

H^'^'ix) = 0. (38) 

It is quite evident that H^' has a smooth zero mass limit, i.e. there is no collinear 
singularity in the (x, w)-dependent structure function H^'{x,w) in Eq. (PT]). This 
must be contrasted with the corresponding structure function describing the linear 
polarization in the event plane which possesses a collinear singularity (see Eq. ([I^)). 
We therefore anticipate that the linear polarization in the gluon-beam plane will be 
generally smaller than the linear polarization in the event plane and that it tends 
to zero in the zero mass limit. 

We are now in the position to calculate the linear polarization of the gluon in 
the gluon-beam plane. The general expression for the linear polarization in the 
gluon-beam plane reads 

9i,{Hh + Hi) + g,,{Hh + HI) ^^^^ 



15 



where we have retained the H'^' term in the numerator even though its contribution 
is zero. 

In Fig. 8 we show the (x/xmaa;)-dependence of the hnear polarization of the gluon 
in the gluon-beam plane for the top quark case for three different cm. energies. The 
linear polarization is generally small and peaks toward the hard end of the spectrum. 
Also the linear polarization tends to become larger when the cm. energy is increased. 
This is counterintuitive from what was said before about the absence of a coUinear 
singularity for the polarized structure function when f — > 1. The approach to the 
limiting behaviour of P^' (for f — > 1) is, however, so slow that one is still very far 
away from the limiting behaviour for the cm. energies considered in Fig. 8. Note 
that the linear polarization neither vanishes at the soft end of the spectrum nor at 
the hard end. In fact, one obtains 

P^'ir = 0)= 9ni-2i3-v')v + {l-v'){3 + v')U) 

2((3 + v^)gn + 3(1 - v^)g^,) {2v - (1 + v^)U) ^ 

at the soft gluon point, where 



1 - e, u = \n (i±^) . (41) 

In the zero mass limit one obtains P^'(O) ^ as expected. At the hard end of the 
spectrum for x = x^aax = 1 — ^5 the linear polarization is given by 

(3 - v^)gn - (1 -v-^jgu 

Note that the mass zero limit of Eq. (^^ is finite, i.e. -P^'(a;max) — > 1/2 as f — 1. 
Again one would naively expect P^' — ^ in this limit. However, one is picking up a 
phase space zero in the denominator of the polarization expression at the hard gluon 
point which cancels the coUinear singularity. The overall effect of the cancellation 
is that the polarization tends to a finite value as f — > 0. In fact, when one plots the 
energy dependence of the linear polarization for very small quark masses, the linear 
polarization tends to zero over the whole energy range, but shoots up to the value 
0.5 very close to the endpoint. 



6 Summary and conclusions 

We have provided a detailed discussion of the linear and circular gluon polarization 
of gluons produced in association with heavy and light quark pairs in e~^e~ anni- 
hilations. We have studied beam polarization and polar orientation effects on the 
polarization of the gluon. The linear polarization of the gluon remains close to its 
classical soft gluon value of 100% and is only mildly dependent on the polarization 
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of the beam, on the polar orientation of the gluon and on the flavour of the heavy 
quark pair produced in association with the gluon. Quark mass effects enhance the 
linear polarization of the gluon when compared at the same scaled enery scaled to 
the maximal energy. The circular polarization of the gluon is strongly dependent on 
the polarization of the beam, on the polar orientation of the gluon and on the flavour 
of the heavy quark pair produced in association with the gluon. The dependence 
on these effects can be exploited to optimally tune the circular polarization of the 
gluon. However, to see a nonzero circular polarization of the gluon one needs to tag 
the flavour of the associated quark/ ant iquark since the circular polarization is odd 
under charge conjugation and averages to zero without flavour tagging. This should 
not be so difficult for the heavy-flavoured quarks. To actually measure the circular 
polarization of gluons one needs to study the fragmentation of circularily polarized 
gluons into polarized particles whose polarization can be measured. 

If one aims to study gluon polarization effects in the splitting process e^e~ — >■ 
tiG{G GG,qq), the present on-shell calculation should be sufficient to identify 
and discuss the leading effects of gluon polarization without having to perform a 
full 0{a1) calculation of e'^e" —>■ ttGG and e~^e~ ttqq 
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Figure Captions 



Fig. 1: Quark-gluon opening angle distribution for soft gluons illustrating the 
"dead-cone" effect 

(a) opening angle distribntion at a fixed cm. energy of — 500 GeV 
for gluon energies of 0.5, 1, 2, and 3 GeV 

(b) opening angle distribution for a fixed gluon energy of 6 GeV at 
three different cm. energies of 370, 500, and 1000 GeV 

Fig. 2: (a) energy dependence of the linear polarization of the gluon in the 
event plane 

(b) energy dependence of the differential cross section e^e" ttG 
and e+e~ — >■ ccG 

Fig. 3: (a) energy dependence of the linear polarization of the gluon in the 
event plane for different cm. energies 

(b) energy dependence of the differential cross section for e+e^ ttG 
for different cm. energies 

Fig. 4: (a) energy dependence of the circular polarization of the gluon pro- 
duced with a top quark pair and a zero mass up-type quark pair for 

different cm. energies 

(b) energy dependence of the circular polarization of the gluon pro- 
duced with zero mass up-type and down-type quarks and with a bot- 
tom quark pair at = 1000 GeV 

Fig. 5: (a) energy dependence of the circular polarization of the gluon in 
e'^e~{hJ) — > ttG with electron beam polarization 
(b) energy dependence of the differential cross section e'^e~{hJ) — > ttG 
with electron beam polarization 

Fig. 6: (a) polar angle dependence of the linear polarization of the gluon in 

the event plane 

(b) polar angle dependence of the circluar polarization of the gluon 

(c) polar angle dependence of the differential cross section 

Fig. 7: (a) Mean hnear polarization of the gluon in the event plane plotted 
against the cm. energy for different cutoff energies 
(b) total cross section for e^e~ ttG plotted against the cm. energy 
for different cutoff energies 

Fig. 8: Energy dependence of the linear polarization of the gluon in the gluon- 
beam plane 
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Figure 1(b) 
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Figure 2(a) 
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Figure 2(b) 
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Figure 3(a) 
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